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SUMMARY 


The paper summarizes results of an investigation of surface pressure fluctuations and response of panels 
underlying attached and separated turbulent boundary layers and shock waves at NASA Ames Research Center. Exten- 
sive tests of a large assortment of axisymmetric and two-dimensional models have been conducted at transonic and 
supersonic Mach numbers to 3.6 to study the pressure fields. Assorted fixed-edge flat panels have been tested at 
Mach numbers from 1. 6 to 3. 6 in attached and completely separated flow fields and also in mixed flow with a step 
induced shock wave oscillating on the panels. The surface pressure fluctuations are described in terms of broad- 
band rms, spectral density, and spatial correlation information. The effectiveness of parameters for scaling the 
pressure fluctuations is also illustrated. Measurements of the amplitude and strain response of the panels are com- 
pared with response computations by the normal mode method of analysis. 


NOTATION (Note: All dimensions in standard metric units) 


A 

surface area of the panel 

h 

panel length in x direction 

A 

m 

normalization factor for $ 

m 

A 2 

panel length in x^ direction 

B~ Eh3 
12(1 - V 2 ) 

bending stiffness 

m 

Of 

generalized mass for the fyth mode 



M =o 

free stream Mach number 

c 

damping coefficient 

P 

local static pressure 

C 

static pressure coefficient, (p - pj/q^ 



P 


p(x,t) 

pressure fluctuation acting on the panel 

d 

thickness of panel 


free stream dynamic pressure 

D 

diameter 

q (t) 
a 

generalized coordinates 

E 

Young's modulus 




Q „<«> 

Fourier transform of q 

f 

frequency 

Or 

Of 



S(x\x", to) 

two-sided cross spectral density function 

f ft) 

generalized force 



Of 


S(x, to) 

two-sided power spectral density function 

fjt/U 

Strouhal number 



c 


s 

flow separation length measured upstream 

F, * 4U/V 
1 1 c 


s 

from x g = 0 

=» 4f l/U 
2 2 c 


t 

time 

F 

Of 

Fourier transform of f 

OL 

u 

velocity 

G(t 4 =* 2S(u>) 

one-sided power spectral density function 

u c (f) 

narrow-band convection velocity 

G(C, co) = 2S& co) 

one-sided cross spectral density 

W(K, t) 

displacement normal to the panel surface 

h 

shoulder height of cone frustum or two- 

W(x, ui) 

Fourier transform of w 


dimensional wedge 

X 

s 

1 ong itud in al di s t an c e me asur e d ups t re am 

H (co) 

frequency response function 

from shoulder of cone frustum or two- 

Of 



dimensional wedge 

j <<*>> 

real part of longitudinal acceptance 



mr 


X 1 

location on the panel in the longitudinal 

j f <^) 
ns' 

transversal acceptance 

(or streamwise) direction 

k <W) 

negative of imaginary part of 

x 2 

location on the panel in the lateral 

mr 

longitudinal acceptance 


direction 
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. X J 

coordinate referring to location on 

e 

angle of cone frustum or wedge 

L Z 

y i = vT 

the panel 

p 

correl ation coefficient 


* a (£) 

mode shape function of the panel 

y 2 =X 2 Ae 2 



mode shape function, of a panel with unit 

xAyy 2 ) 

normalized coordinates referring to 
location on the panel 

dimensions 




CO 

angular frequency 

6 

boundary layer thickness 

to 

Of. 

natural frequency of the oth mode 

6* 

displacement boundary layer thickness 




0 

phase angle of cross spectral density 


Krone eke r delta 


function 


ry jS 


jL.,±± t jL 

.4 <2.2 ^ 4 

oXj Bx 1 3 x 2 dx 2 


Subscripts 

norm normalized quantity 


| = x 1 - x” coordinate referring to the separation 
distance of two points on the panel 


Of =» (m,n) mode index 
3 = (r, s) mode index 


^2 “ y 2 ' y 2 

ncn v v 2 ) =x’ - 1” 

€ strain 

/i mass per unit area of the panel 

p Poisson ratio 

v loss factor for the ctth mode 

a 

Or attenuation coefficient 

y coherence function 


o quantity evaluated immediately ahead of 

detached shock wave 

00 quantity evaluated at the free stream 

1 streamwise direction 

2 lateral direction 

p quantities related to the excitation field 

d quantities related to pane] displacement 

c quantities related to strain 


1. INTRODUCTION 

The pressure fluctuations in regions of attached and separated turbulent boundary layers and shock waves 
adjacent to surface of aerospace vehicles give rise to structural vibrations throughout atmospheric flight. The study 
of these excitations and resulting vibrations is of importance in determining stress, fatigue life of structures, and 
noise transmission into the interior of the vehicle. Unfortunately, the analysis of this type of vibration is complicated 
by the inherent random characteristics of the excitation pressure fields, and the difficulty of analytically describing 
the vibration of a realistic structure. For these reasons, early investigators of this problem considered only a 
hypothetical flow field and made the simplifying assumption that the structure, almost invariably either a beam or a 
rectangular panel, is infinitely large (refs. 1-3). This assumption gives rise to a solution in terms of the mean square 
displacement of the panel as a whole, but not the displacement as a function of location on the panel. 

Investigations dealing with the excitation fields for the attached and separated turbulent boundary layers 
and regions of mixed flow, including oscillating shockwaves, have been numerous (the references of l'efs. 4-6 yield 
an appropriate lengthy bibliography). A review of the literature indicates, however, that the required statistical 
information to describe the surface-pressure fluctuations is reasonably complete only for subsonic attached turbulent 
boundary layers. Investigations at transonic and supersonic speeds have primarily included measurements of pressure 
fluctuation intensities and very limited analysis of power spectra and/or spatial correlation. With the exception of 
results in references 4-7, these latter forms of analyses have only been published for attached flow. 

With respect to response, recent investigations have considered finite-size rectangular panels underlying 
subsonic attached turbulent boundary layers, but a simplifying assumption was made that the panels were simply sup- 
ported (refs. 8-10). This assumption simplifies the algebra tremendously, permitting solutions to be expressed in 
closed forms; and although these analytical results give better agreement with experiment than those obtained with 
the infinite panel assumption, they tend to overestimate the response of a realistic panel. Previous analyses also 
tailed to give an accurate prediction of the wave matching between the flexural wave of a panel and the pressure wave. 
The estimation of this wave matching is important particularly when the matching occurs at one of the resonant fre- 
quencies of the structure and this causes a large structural response. With the advent of modern high speed computers 
a theoretical analysis of a finite rectangular panel with clamped edges under the excitation of a turbulent boundary 



5-3 


layer is now feasible. Analytical integration in closed forms are not necessary, since numerical integration can be 
carried out with no algebraic simplification of the integrand. This digital -computer -“or tented approach has the follow- 
ing additional advantages over an analytical approach: (1) The transparency of the problem is preserved* as very 
often the physics of the problem is lost amongst a great length of closed-form mathematical formulae; (2) future 
developments are simplified as the basic computer program can be modified to describe different flow fields and/or 
different structures. 

With the deficiencies in both the input and response problems in mind, a research program has been under- 
taken at NASA, Ames Research Center, to improve the statistical description of the random pressure fluctuations 
underlying attached and separated turbulent boundary layers and shock waves; and to improve the analytical capability 
for computation of the displacement and strain of realistic panel structures when excited by each of these flow fields. 
The investigations of pressure fluctuation inputs and response have spanned the subsonic, transonic, and supersonic 
speed ranges up to a Mach number of 3.5. 

Sketches of models used for the investigation of pres sure -fluctuation inputs are shown in Fig. 1. The 
basic configurations were 0. 0508 m and 0. 254 m diameter ogive cylinders and the Ames 9- by 7-foot and 8- by 7-foot 
supersonic wind tunnel (SWT) walls to investigate attached turbulent boundary-layers. Cone frustums ahead of axi- 
symmetric rings and two-dimensional wedges of different heights and a variety of angles from 15° to 90° were added 
to the cylinders and walls to investigate regions of separated flow and shock waves. The different model sizes and 
tunnel walls provided large variations of the thicknesses of both the attached and separated boundary layers to investi- 
gate scaling relationships to establish the most effective parameters for nondimensionalization of the pressure fluctua- 
tions. The thickness of attached boundary layers, for example, varied from approximately 0. 00406 m on the 0. 0508 m 
diameter model to 0. 135 m on the wind-tunnel wall. Separated -flow lengths ahead of 45° frustums and a 0. 2032 m 
high 45° wedge varied from 0. 041 m to 0. 89 m. 


ATTACHED TURBULENT 
BOUNDARY LAYER 


.0508 m diam 


.254 m diam 



2- DIM {W.T. WALL) 


SEPARATED FLOW AHEAD 
OF RAMPS AND STEPS 



Figure 1. Configuration investigated for study of surf ace -pres sure fluctuations. 

The surface pressure fluctuations were measured with two different sizes of semiconductor type pressure 
transducers having sensing diameters of 0. 00279 m and 0. 00121 m. The transducers were mounted flush or slightly 
submerged to the surface in longitudinal, lateral, and diagonal arrays. In addition, detailed static-pressure distribu- 
tions and boundary-! aye r measurements were obtained. More details on the test procedure, data acquisition and 
analysis, and overall scope of the pres sure -fluctuation investigation can be found in references 4 and 7. 

For the response part of the research, the first step in the program was to develop the analytical capa- 
bility and a practical computer program for computation of the amplitude response of a clamped-edge panel. The 
excitation of a subsonic attached turbulent boundary layer was considered first because of the availability of corres- 
ponding excitation data by Bull (ref. 11), and response measurements by Wilby (ref. 10). The method of analysis 
(method of normal mode used together with the technique of spectral analysis) and comparisons with Wilby 1 s response 
data are shown in reference 12. These initial results show favorable comparisons between the computed and measured 
response. 


The analysis and computer program have subsequently been expanded to include the amplitude and strain 
response of simply-supported and clamped-edge panels excited by supersonic attached and separated turbulent boundary 
layers. Corresponding experiments have been conducted to measure the amplitude and strain response of a variety of 
panels with different length -to -width ratios and thicknesses. Panel lengths and widths varied from 0. 1524 m to 
0. 3045 m and thicknesses varied from 0, 5588* iG ~ 3 m to 2. 234X1CT 3 m. The tests were conducted in the same facilities 
(9- by 7 -foot SWT and 8- by 7-foot SWT) used for. the supersonic pres sure -fluctuation studies between Mach numbers 
of 1. 6 and 3.5. An illustration of the test set up for separated flow is shown in Fig, 2. The panel response tests were 
conducted with the 0. 2032-m high 90° step, whereas pressure fluctuations were also measured ahead of 23° and 45° 
wedges. Separated flow lengths ahead of the steps varied from approximately 0.31 m to 1. 10 m depending on Mach 
number and the attached boundary layer thickness. 

Three test conditions with the panels underlying attached flow, separated flow, and mixed flow with an 
oscillating shock wave were investigated. Attached flow was obtained by removing the step. Separated flow was 
obtained with the step positioned as shown in Fig. 2 so that the panel was about one-third the distance between the 
shock wave and the step. The step was moved further downstream to position the shock wave on the panels. 



5-4 



Figure 2. Panel -response test installation. 
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Panel displacements were measured with non- 
contacting capacitance probes that could be remotely posi- 
tioned to different panel coordinate positions. A steady -state 
cavity pressure was carefully maintained at wall static 
pressure at the streamwise midpoint near the lower edge of 
the panels for all test conditions. 

The analyses of all the results of the pressure 
fluctuation and panel response investigations are not com- 
plete. However, a large amount of data have been amassed 
for attached and separated flows at supersonic speeds, and 
therefore these data and the corresponding panel response 
measurements and analyses are the subject of this paper. 

2. PRESSURE -FLUCTUATION INPUTS 

2. 1 Broadband Characteristics 


The surface-pressure fluctuations measured in 
attached flow on the walls of Ames wind tunnels are shown 
in Fig. 3 along with comparative data from a few other 
investigations . The data presented are considered by the 
authors to most reliably represent the broadband intensities 
of the pressure fluctuations, /p2/qco, for the Mach number 
range from 0 to 3.5. Many more data points could be added 
to the figure from other earlier investigations, but the 
spread of data would then be increased markedly. Most of 
the other available data considered have been rejected 
because of uncertain effects of transducer size, installation, 
and extraneous noise. The data from the Ames 9- by 7 -foot 
SWT and 8- by 7 -foot SWT and the XB-70 represent measure- 
ments at several locations on a rigid plate in the tunnel wall 
and on the aircraft. The spread in data of about 30% (3 dB) 
is typical and is believed due to a nonuniformity in the 
unsteady pressure field. The fact that data can be con- 
sistently repeated within 0. 3 dB when different transducers 
of the same kind are installed at a specific location on the 
wind tunnel wall supports this argument. Generally, the 
wind tunnel and flight data in Fig. 3 co mpa re favorably, 
and both show about a 50% decrease in'/p^/q O0 between 
Mach numbers of 1 and 2.5. At M > 2. 5 the available wind 
tunnel data show a tendency for J p^/q ro to increase slightly 
from a minimum of about 0. 0035 at M = 2. 5 to about 0. 004 
at M = 3.5. 
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Figure 3. 


Broadband characteristics of pressure 
fluctuations underlying attached 
turbulent boundary layers. 


separation points were not scaled by the x Q /h parameter. 


Typical longitudinal distributions of mean static 
and fluctuating pressures illustrative of the separated flow 
on the models of this investigation are shown in Fig. 4 (the 
effect of diameter and frustum height in Fig. 4(a) and the 
effect of Mach number on the flow field in the vicinity of 
the response-test panel is shown in Fig. 4(b)). As pre- 
viously described (refs. 4 and 15), the mean static pressure, 
initially at free -stream static pressure, increases rapidly 
as the flow encounters the detached frustum shock near the 
separation point and then approaches a plateau as the transi- 
tion from attached to separated flow is completed. The 
static -pres sure coefficients in the region of the plateau 
were nearly the same on all the axisymmetric models, but 
were approximately 20% lower than the corresponding meas- 
urements on the two-dimensional wall. The locations of 
the static-pressure rise associated with the shocks and 
Further study of the shock separation distances has shown 


that the differences are due to the variation of h/6 Q on these models rather than to Reynolds number effects. 


The characteristic feature of the broadband pressure fluctuations (Fig. 4(a)) can be related to features of 
the mean static-pressure distributions. The pressure fluctuations increase rapidly at the shock, as identified by the 
rapid pressure rise, and reach a maximum intensity where the slope of the static-pressure curve is maximum. In the 
fully separated region, the fluctuating pressure descends from this maximum to a plateau level. Correspondingly in 
this same region the static pressure descreases toward a plateau level. It can be noted that the plateau level was 
relatively independent of the shoulder height of the frustum that caused the separation on the 0. 254-m-diameter models. 
The lower fluctuating -pres sure intensities on the 0. 0508-m -diameter model were found to be due to an insufficient fre- 
quency range of the recorded data. The higher fluctuating pressures on the two-dimensional model are associated 
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a) Effect of diameter and frustum/wedge height. 


b) Effect of Mach number. 


Figure 4. Longitudinal distribution of steady and fluctuating static pressures in region of separated flow. 

with the higher static pressures in the separated flow region on the two-dimensional model. Generally, the Mach num- 
ber effects on the broadband pressure fluctuations (Fig. 4(13)) in the separated flow are not significant, however, there 
is a trend showing a slight increase in intensity with increasing Mach number. 


2.2 Power Spectral Densities 

Power spectra representative of pressure fluctuations in the vicinity of supersonic separated flow are 
shown in Fig. 5. The curves* show the distinct differences in shapes and mean-square amplitudes that distinguish 

separated flow spectra from spectra measured in regions 
of the nearby shock wave and attached turbulent boundary 
layer. Shock spectra are characterized by a very steep 
slope and high intensities at low frequencies. The lowest 
frequencies in these data are 10 Hz. The separated-flow 
spectra were at least one decade lower than the shock 
spectra at low frequencies. Having a less steep slope, 
however, the separated-flow spectra crossed the shock 
spectra, and were the highest spectral intensities obtained 
at intermediate and high frequencies. The power spectra 
of the pressure fluctuations beneath the attached boundary 
layer were more than two decades lower than the separated- 
flow spectra at low frequencies. Although all the spectra 
appear to be converging toward a common level and slope 
at higher frequencies, there was still at least one decade 
difference between the attached and separated flow spectra 
at the highest frequencies of this investigation. The spectra 
labeled 1 ’intermittent shock” illustrate the effect of an 
intermittent excursion of the shock to the maximum limit 
of its upstream oscillation. Only the very low frequency 
components of the pressure fluctuations are increased in 
intensity, while the intermediate and high frequencies retain 
the intensities of the attached boundary layer. 

2 , 2 . 1 No ndimensionalization of Power Spectra 

Many geometric and flow parameters were 
considered for the nondimensionalization of pressure fluc- 
tuations when presented in the frequency domain. Various 
combinations of length and velocity were applied to the 
generally accepted reduced frequency parameter having 
the form (frequency x length) /velocity. The mean- 
square amplitudes per unit reduced frequency were 
nondimensionalized by (pressure) Some of the 
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Figure 5. 
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Typical power spectra of pressure fluctua- 
tions underlying supersonic flow. 
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reduced -frequency parameters investigated were: fS /U o , f6*/U roS f S /U, f(S^ - x )/U , fd/U » and 
f6/U « 


It became evident early in the investigation that there was no reduced-frequency or spectral “dens tty 
parameter that was the most effective for scaling the data from the many models,, particularly for separated flow, 
for the full range of frequencies considered, (Certain characteristics of the pressure fluctuations were found to be 
different within different ranges of reduced frequencies.) It was clear however, that the parameters ffi/U and 
f(S g - x s )U involving local bound ary “layer thickness or flow-separation dimensions and local velocity were best for 
scaling frequencies in separated flow. In this paper, since both attached and separated boundary layers are considered, 
the nondimensionalization of spectra are illustrated in terms of f6/U. The data used for the response computation are 
also presented in terms of f6/U« since the local velocity in the separated region is frequently not available and also 
since differences in scaling with U vs Uco are relatively minor. Also minor differences were found in the effective- 
ness of free -stream versus local dynamic pressure for scaling the spectral density; therefore, the choice has been 
made in favor of the more conveniently available qoo. 


Local boundary layer thicknesses used for scaling the pres sure -fluctuation data were based on measure- 
ments of boundary-layer profiles for all the attached-flow test cases and on profiles measured within the separated 
flow region on the d - 0 . 254 m, h = 0 . 0508 m, ft = 45 ° model. The separated-flow boundary layer measurements are 
shown in Fig. 6 to illustrate that the boundary-layer growth rearward from the separation point was sufficiently linear 

to allow a simple linear interpolation of local boundary- 


0=45° h = . 0508 m 



layer thickness. Local boundary-layer thicknesses in 
separated flow were therefore estimated for all the models 
from the relationship 6 = 6 Q + h (S g 


■ x )/S 
s s 


Figure 6. Boundary-layer thickness 
of separated flow. 


Typical scaling of spectral measurements 
obtained on the different models at M = 2. 0 is illustrated 
in Fig. 7. The spectrum selected for each model repre- 
sents an approximate mean of from 10 to 20 measurements, 
and in most of the separated flow cases the mean spectrum 
was obtained near the center of the separated flow region. 
The closeness of fit of the spectra scaled by the reduced 
frequency parameter f6/U (Fig. 7) illustrates about the 
best attainable collapse of data. Most of the separated 
flow results fit within a factor-of-two spread on the mean- 


squared spectral density scale. This is considered a good fit of random dynamic data obtained on so many different 
models, and indicates relative independence of the pressure-fluctuation characteristics on geometry. 


ATTACHED FLOW 


SEPARATED FLOW 



M<o=2.0 
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Figure 7, Scaling of power spectra. 


2.2.2 Mach Number Effects of Power Spectra . 

The variation of power spectra in attached and separated flows between Mach numbers of 1. 6 and 3. 5 is 
shown in Fig. 8. Each of these spectra represents an average of from 10 to 20 spectra measured on a rigid plate 
installed in the "panel-response” test fixture in the walls of the Ames 9- by 7-foot SWT and 8- by 7-foot SWT, These 
data, form the basis of the representation of the excitation spectrum used for panel response calculations. Note that 
the data are nondimensionalized by free-stream velocity. 

The power spectra for attached flow generally are very similar for all Mach numbers with the exception of 
a slightly lower measured spectrum at M = 3. There is a reasonably consistent trend at low reduced frequencies, 
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Figure 8. Variation of power spectra with Mach number. 


f6/U OT < 0.05, that indicates a reduction in spectrum level with increasing Mach number. These variations should be 
taken lightly, however, because of reduced statistical accuracy of analysis at the lower frequencies. In contrast with 
the attached flow ease the power spectra measured in the separated flow show a very consistent variation with Mach 
number. The spectral levels decrease with increasing Mach number at f6/U m <0.3 and increase with increasing 
Mach number at f6/U w >0.3. These results illustrate the separation of frequency regimes at fd/Ueo = 0. 3. In each 
regime the pressure fluctuations are dependent upon different phenomena. The surface pressure measurements in the 
lower frequency region are mainly sensitive to the transmitted fluid fluctuations in the turbulent flow of the separated 
boundary layer; the radiated energy is negligible. In the higher frequency region, it is believed that the surface pres- 
sure fluctuations are strongly influenced by an eddy Mach wave radiation phenomena. This phenomena was investigated 
by Ff owes -Williams (ref. 16), and although his results do not apply to the study of surface pressure fluctuations a study 
currently in progress by Dr. Richard D. Rechtxen of the University of Missouri, Rolla, Missouri, (NASA-Ames Grant) 
indicates a dependence of these data on Mach wave radiation. 


2 . 3 Spatial Correlation of Pressure Fluctuations 


Chyu and Hanly (ref. 4) presented some data 
from this investigation that showed the power and cross 
spectra and space -time correlations of fluctuating pres- 
sures underlying the supersonic attached boundary layer 
on the 0,254-m diameter ogive cylinder and the separated 
boundary layer ahead of the 0. 0508-m high frustum. 

Among other things, they investigated the co- and quad- 
spectral density and coherence as functions of a wave num- 
ber parameter f£^/U c (f). Illustrative coherence functions 
are shown in Fig. 9. The r*e suits indicate a similarity 
between the coherence measurements in attached and 
separated flow at supersonic Mach numbers, but only if 
the velocity term is U c (f). The use of convection velocity 
instead of free -stream velocity is necessary, since for 
separated flow, U c (i) varies significantly with frequency. 
The general trends of the data are also similar to the 
subsonic-flow measurements by Bull (ref. 11). It can be 
seen that the envelope of the coherence functions for vari- 
ou s sp at i al d i s t ance s 1 1 / 6* decreases exp one n t i ally wi th 
increasing frequency, however, the envelope unfortunately 
represents the coherence only at high frequencies. 

Chyu and Au-Yang (ref. 12) used BuIUs 
coherence data to represent the subsonic attached- 
boundary-layer spatial correlation for application to the 
response computation. The fact that the coherence func- 
tions do not collapse complicates the empirical expression 
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Figure 9. Longitudinal coherence of pressure 
fluctuations. 
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of correlation, and therefore, a different approach has more recently been investigated. As suggested by Rechtien 
the normalized moduli of the cross spectral densities, jG^£, f)j norm ,^ for available or selected transducer s pacings 
have been curve -fitted to the exponential function |Gp<i£, w| norm “ by the method of least squares to obtain a non- 
dimensional attenuation -Co efficient function <y{4, f)6 0 Typical moduli of the cross-spectral densities of the pressure 
fluctuations in separated flow that were used for the evaluation of the attenuation coefficient are shown in Fig, 10. 


Figure 11 shows the attenuation-coefficient function measured in separated flow on the 0. 254-m diameter 
model with the 0. 0508-m high 45° frustum shoulder to illustrate its effectiveness of describing spatial correlation. 

The function is shown to be reasonably independent of the number of moduli, the transducer spacing, or the reference 
locations involved in the curve fitting. These results indicate that the flow is relatively homogeneous within the limited 
area of the separated flow region where the measurements were obtained. It appears that the attenuation coefficients 
decreased slightly at higher frequencies as the larger transducer spacings were used in the analysis. The variation is 
considered insignificant, however, in light of the state of the art of the structural part of the problem of predicting 
response to random turbulence. Attenuation coefficients obtained from transducer arrays oriented longitudinally, 
diagonally, and laterally to the free-stream flow indicate that the decay of correlated turbulence was independent of 
orientation at f6/U < 0. 06. The predominant turbulence is therefore nonconvective at the lower frequencies, and con- 
tours of equal spatial correlation would be circular. At f6/U > 0. 06 the attenuation coefficients were progressively 
lower as the angularity of the transducer orientation changed from lateral to longitudinal, indicating extended correla- 
tion in the direction of free-stream flow. 
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Figure 10. Variation of normalized cross spectra 
with longitudinal spacing of transducers 
in separated flow. 
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Figure 11. Attenuation coefficients measured on the 
0. 254-m diameter model in the region 
of separated flow, h = 0. 0508 m, 

0 = 45°. 


A comparison of longitudinal attenuation coefficient functions obtained on the different models in separated 
flow is shown in Fig. 12 to illustrate the effectiveness of a 6 for scaling for separated flow. The results indicate that 
the nondimensionalized attenuation coefficients were in relatively good agreement, showing a total spread between all 
the curves of generally less than a factor of two. 

Typical longitudinal attenuation coefficients at Mach numbers from 1. 6 to 3.5 are shown in. Fig. 13 for both 
attached and separated, flow on the wind tunnel walls. As with the spectra (Fig. 8) the separated flow data show a more 
consistent variation with Mach number than the attached flow data, particularly at the lower frequencies. The results 
show that the pressure fluctuations in attached flow were most highly correlated at r6/U ro « 0.4; whereas, in separated 
flow the maximum correlation occurred at the lowest frequencies. The dips in separated flow o^6 curves in the fre- 
quency range 0. 3 < f6/U ro <1.3 indicate increased correlation that is believed to result from the Mach wave radiation 
phenomena discussed in connection with the observed Mach number effects on power spectra. 
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Figure 12. Scaling of attenuation coefficients. 


The convection, of turbulence in attached and 
separated boundary layers results in an angle of phase 
between spatially correlated components of pressure. Since 
this- phase angle can result in matching of panel flexural 
waves and pressure waves and thus cause an enhancement 
of response, it is necessary that it be accounted for in 
response computations. Convection of the surface -pressure 
fluctuations can be studied directly in terms of convection 
velocities or in terms of the phase angles of the cross 
spectra. The broadband convection velocity derived from 
space -time correlations is the most convenient description 
of convection characteristics for attached flow, but Chyu 
and Hanly (ref. 4) showed that it is inappropriate for 
separated flow since convection velocities vary with 
frequency. 
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Figure 13. Effect of Mach number on longitudinal attenuation coefficients. 

Typical phase angles that result from convection of the turbulence in separated flow are shown in Fig. 14. 
The data were obtained on the 0. 254-m diameter model with axi symmetric flow and show the effects of the angular 
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Figure 14. Typical phase angles of cross spectra of 
pressure fluctuations underlying 
separated flow. 


orientation of a transducer array to the free -stream flow 
direction. At low frequencies, f < 300 Hz, the phase angles 
were zero for all angular orientations. A zero phase angle 
indicates a zero or infinite speed of convection of the pres- 
sure fluctuations between two points, as would exist with a 
stationary pulsating bubble or for normal incidence of an 
acoustic radiation process. The equivalent of the pulsating 
bubble could occur if the turbulence is related to fore-and- 
aft oscillations of the detached shock wave. The phase- 
angle measurements generally indicate an absence of lateral 
convection in the separated flow over the full range of fre- 
quencies investigated. This result was evidenced by the 
near zero lateral measurements and by the agreement 
between longitudinal and diagonal measurements of 9/^ 
(phase angle per meter of longitudinal component of trans- 
ducer spacing). 

2 * 5 Narrow-Band Co n vection Velocities 

When significant phase angles are measured, 

they can be converted to convection velocities, as defined 
by 360 f(£/0), and thus reveal the narrow-band convec- 
tion characteristics of the turbulence. Figure 15 shows that 
on the 0. 254-m diameter model, the convection velocities 
of the pressure field in separated flow varied from a mini- 
mum of approximately 0. 2 U at f6/U ^0.6 to a plateau 
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Figure 15. Typical narrow-band convection velocities. 


near 0. 8 U at £6/U 0. 8. The variation of convection, 

speed with frequency infers that the predominant turbulence 
at different frequencies between. 0. 03 < f6/U < 0. 8 is gen- 
erated at different levels of the boundary layer ranging from 
slightly above the zero velocity line to the free-shear layer. 
Upstream convection was not detected between points within 
the separated flow region. The turbulent eddies generated 
in the reverse flow region apparently contribute little to 
the wall-pressure fluctuations. With the exception of the 
low frequency region of f6/U < 0. 06, the convection 
velocities were relatively unaffected by transducer spacing 
or reference location. It is also shown that they can be 
effectively scaled by the reduced frequency parameter f6/U. 

Representative narrow-band convection velocities 
in attached and separated flows on the walls of the 9- by 7 -foot 
SWT and 8- by 7 -foot SWT at the location of the panel response 
test fixture are shown in Fig. 16. It can be seen that for 
attached flow there was no significant variation in the ratio 
of convection velocity to free-stream velocity 
(U c (f)/U« » 0.75) with frequency or Mach number. The 
variation of convection velocity with frequency in separated 
flow was previously discussed in connection with Fig. 15. 

This variation with frequency was relatively unaffected by 
Mach number at f6/U ro < 0. 6. The plateau values of 
nearly constant U c /Uce at f6/U txl > 1. 0, however, varied 


inversely with Mach number (with the exception of the M OT = 3.5 data) so that the absolute convection velocity, U c , 
was nearly constant. The convection velocity in this plateau region varied from approximately 354 m/sec (1150 fps) 
at Mo, = 1. 7 to 366 m/sec (1200 fps) at Ma> = 3. 0. The reduced frequencies where the beginning of the plateau occurs 
can be identified with changes in the characteristics of the power spectra and attenuation coefficients. 
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Figure 16. Variation of narrow-band convection velocities with Mach number. 

2. 6 Representation of Surface-Pressure Fluctuations 

The fundamental information needed to describe the excitation of a structure underlying a turbulent flow is 
the cross spectral density of the fluctuating pressures. 


In a homogeneous turbulence the cross spectra can be expressed in the form 




= P~<yL,,ye 


- iW V U c 


( 1 ) 


where °°)| * s cross correlation coefficient between two points separated by distance 4^ and ^ * n the 

Xj and x| directions. Experimental data have also shown that the correlation coefficient can be separated in coor- 
dinates as follows: 




(2) 



5-11 


where |p~(? lS 0, to)| = exp [- 0 ^] and jo_(0, a))| = expC-a^] # 

The convection velocity U c is related to the phase angle 0 of the cross spectral density of the fluctuating pressure by 
U c * 360 f^/0. For this investigation, It has been chosen to represent U c in terms of 0 since U c varies significantly 
with f for separated flow and since 0 -* 0°, U c - «. it should also be noted that 0/^ is an independent function of f 
and is therefore not scaled by a reduced frequency parameter. 

To derive expressions that represent the selected functions used to describe the fluctuating pressures in 
supersonic attached and separated turbulent boundary layers, the large amount of experimental data previously dis- 
cussed have been averaged at each test Mach number (1. 6, 2.0, 2.5, 3.0, 3.5) and approximate means of the envelopes 
of these data have been curve-fitted as shown in Fig. 17. 
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Representation of surface-pressure fluctuations. 
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The empirical formulae that express the power spectrum, attenuation coefficients, and phase angles for 
attached and separated flows are as follows: 


a) For attached flow 
Power Spectrum 


G_U /q 2 6 = exp f-12. 470 - 0. 639x - 0. 263x 2 4 0. 015x 3 4 0, G17x 4 4 0. 002x 5 ] 

p co *a 3 1 j 


( 3 ) 


Attenuation Coefficients 


a 6 “ exp {0. 060 + 0. 941x 4 0. 380x 2 - 0. 023x 3 - 0. 020x 4 - 0. 002x 5 } 
0^6 = exp {2. 163 + 1. 29 lx + 0. 262x 2 - 0. 132x 3 - 0. 043x 4 - 0. 0033x 5 } 


(4) 

(5) 


Phase Angle 


0/4 = 39. 37 {115. 14 4 131. 99x + 69. 89x 2 + 18. 84x 3 4 2. 44x 4 4 0. 12x 5 } deg/m 


( 6 ) 


b) For separated flow 
Power Spectrum 


G^U/q 2 6 - exp {-8. 094 - 1. 239x - 0. 295x 2 - 0. 090x 3 - 0. 014x 4 - 0. OOlx 5 } 


(?) 


Attenuation Coefficients 


a 6 = exp {1. 031 4 0. 666x 4 0. 645x 2 4 0, 270x 3 4 0. 043x 4 4 0. 002x 5 } 
0^5 - exp [1. 797 4 1. 239x 4 0. 536x 2 4 0. 122x 3 4 0. 010x 4 ] 


( 8 ) 

(9) 


Phase Angle 


0/^ = 39. 37 {90. 54 + 89. 28x + 50.47x 2 + 15. 69x 3 + 2. 32x 4 + 0. 13x 5 } 


( 10 ) 


Where for both a) and b) above: x = log^fS/U^) 

3. RESPONSE OF FLAT CLAMPED-EDGE PANELS 
3. 1 Method of Analysis 


The displacement w(x, t) of a vibrating panel 
(Fig. 18), is assumed to obey the classical thin plate 
equation 


/i w 4 cw 4 BV w - p (x, t) 


( 11 ) 


where d and B are constants and independent in the 
present analysis. It is assumed further that w(x, t) can 
be expanded in terms of 0 (x) as follows: 



w(x,t) 


(12) Figure 18. Schematic of a panel with flow excitation. 


Here 0^ is assumed to be properly normalized and separable in coordinates, i.e. , 


* 0 ® ** = 6 a 0 


(13) 


and 


^ a (x) = 


( 14 ) 
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It can be shown that for a panel with clamped edges takes on different forms according to whether m is even 
or odd; 

a) If m is odd, 


t/> (x.) 

m i 


A &, 

m v i 


■ cos y 


L i 1 


4 s K cosh y 


m \ f,. 2/ m mV A. 2, 


i X 


(15) 


where y are the roots of the equation 
m 


tan Hr + tanh 

u & 


and 


K = 


3 in (V 2 ) 
/ 2 ) 


in sinh/y 


b) If m is even, 


(xd =- 


A JIT 

mr i 


sin y 


/l i) 

2 ) 


+ K sinh j - ~ 
m \ & 2 


(16) 


where y are the roots of the equation 
m 


m , , m A 

tan — - tanh — = 0 


By using the orthogonality condition of the mode-shape function, the normalizing factors A ffl can be found to be 


where 


m Z 


2 

Z -y + sin y 4- K (sinh y 4 y ) if m is odd 

m m m m mm 

2 

= y - sin y + K (sinh y - y ) if m is even 

' m m m m m 


The generalized coordinates q a (t) satisfies the Lagrange equation 

m q (t) + C q (t) + K q (t) = f (t) 
aVv 1 ; a a axr ' a 


(17) 


where the generalized mass 


■/*« 


/ t!>„ (x)dx = 


K = moo 

Of Of 


f (t) 

Of 


= yV (X, t)ll> a (x)dx 


and 


K V 
rv fy 


a oo 


here v is the combined loss factor for both structural and viscous damping. The present state of art requires that 
v must be measured experimentally. The Fourier transform of the Lagrange Eq. (17) takes the form 


Q« (aJ) = V W) V W> 


(18) 
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where 


H (to) - 
or 


m f (co 2 - o o 2 ) + ip a) 2 ] 
Of a a 


(19) 


Equations (12) and (18) together give 


W(x,w)= ^H a (W)F^(w)!|) a (x) 


( 20 ) 


The displacement power-spectral density is related to W (x_, co) by 


S d (x, co) = Lim - W T * (x, co) W T (x, co) 
T — *00 


( 21 ) 


where is the truncated Fourier transform of w. Equations (17), (18), (20), and (21) together give the displace- 
ment power spectral density in the form: 

v*’» - v'> i ?,« KM 2 


* v'E 


[g ,aV (F i ) + h ^ k . r (F i )] 


2 2 

m m fl (g + h , 
a jS^cyjS at 


(22) 


where 


g = (co 2 -co 2 ) (to 2 - to 2 ) 4 p l/.w 2 co 2 
ap ' <y p a. p a $ 


h A = p to 2 (co 2 - co 2 ) - p co 2 (co 2 - to 2 ) 
a p a a a a p a 


W F i> 


1 y ,! 

- 2 J / LWI'kfi'" 1 *) 


1TF 1 T! 1 

: — dy ! dy M if m 4 r even 
2 *1 


= 0 if m + r odd 


(23) 


and 



1 y" 
r r 2 

/ 


W = 2 

Jj ww 


>’ V —) 


0 0 



= 0 

if n 4 s odd 




1 y" 
/* /* ^ 


' „ 

V'T/ 

k (FI =2 
mr r 

1 I 

p p( 


O O 



= 0 

if m 4 r even 




dy^dy” if n 4 s even 


TTF^i 

sin — ~ — dy^dy^ if m 4 r odd 


(24) 


(25) 


Here 0 (y) are the normalized mode -shape functions, expressed by Eqs. (15) and (16) with l replaced by unity. The 
physical significance of the structural acceptances J mr is discussed by Ref. 12. The pressure correlation coefficients 
are experimentally shown to be related to the cross spectral density of the fluctuating pressure in a homogeneous 
turbulence and can be expressed by 




S~(o>) 

P 


0~(t ,0,<o) 
P 1 


-ito£ /U 

[ p — (0 , C , to)| e 
I P 2 'i 


( 26 ) 


The strains in the x^ and X 2 directions of a rectangular plate are related to the displacement w by 



_ (l/2)d / 1 
” 1 - V 2 S r 2 



(27) 
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where 


. .2 

1 o w 



- .2 

1 _ ,5 w 

r * ‘ 2 

2 ox 


are the curvatures in x^ and x 2 dii'ections. The strain spectrum can he obtained by 

s (X, W) = lim ~ E* (x, 03 ) E (x, 03 ) (28) 

i T-» T T 


where E. is the truncated Fourier transform of e. in Eq. (27), and can be expressed by 
T 1 

Ej&w) =£ ^(03^(03)$* 

O' 


and 


(29) 


E 2 fe«)=EH Q ,(03)F a (U3)$° 


where 


1 - p2 W - 2 2 / 


,« = nzg).g + „ Uk ) 


*2 212 2 
1 - V \d * 2 


Equations (28-30) together with the expression for i|) in Eqs. (15) and (16) give the strain power spectrum 

S (x,f) = AS (f)2>“(x)] 2 |H^(f)| 2 j mm (F 1 )j; n (F 2 ) 
i a 

E *“(x) $?(x) j' (FJ [g fl j (F ) + h -k (F )] 

i w i w J ns 2* o'jSW' 1 a.8 mr V 

— 




m a m e {s J + h J } 


(30) 


(31) 


3.2 Response Measurements and Computations 

As briefly described in the introduction, tests have been conducted to measure the displacement and strain 
response of several panels of different lengths, widths, and thicknesses in supersonic attached and separated flows 
and in mixed flow with a shock wave positioned on the panels. The response tests of the total input-response investiga- 
tion have only recently been completed and therefore the following results are somewhat preliminary. All the panels 
tested were constructed of magnesium having the following properties: 

10 2 
E - 4. 48X10 N/m 

D = 1780 Kg/m 3 

V = 0. 35 

The first panels selected for analysis had the same dimensions with = 0. 3048 m, = 0.2286 m, = 4 /3, but 

two different thicknesses d - 0. 00118 m and d = 0. 00235 m. 

It was recognized prior to the tests that the installation of strain gages on the panels could affect the stiff- 
ness and symmetry of response and therefore could compromise the comparison between measurements and computa- 
tions. Since strain is a key element in fatigue, however, it was considered more important to include the strain gages 
to verify the strain computations. It would have been desirable to conduct separate displacement response tests, but 
this was economically impractical. 

3.2.1 Systems Damping 

In Eq. (19) v r is the combined loss factor for the complete structural air stream system. The loss factor 
therefore takes into account the effect of hysteretic damping as well as viscous damping. The latter includes the 
interaction of panel motion on the flow field within and outside the boundary layer, and the effect of acoustical radiation 
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into the interior of the structure (cavity) and into the exterior flow field. Although the development (by Dowell, ref. 17) 
of an analytical method for the determination of the aerodynamic damping part of the loss factor is in progress, the 
present state of art requires that V - must be estimated from experience or measured experimentally. 

The loss factors used for the computations of panel response as presented in this paper have been measured 
from the autocorrelations of the response -time histories for the specific panels of interest. Representative loss -factor 
measurements for a 0. 3048 -m long, 0. 2286-m wide by 0. 002 35 -m thick magnesium panel are shown in Fig. 19 for both 

attached and separated flows. The results indicate that only 
the damping of the first mode is significantly higher than the 
damping of all other modes. The first mode damping was 
slightly higher in separated flow than in attached flow, but 
all other modes were relatively unaffected by the flow condi- 
tion. The effects of Mach number are not clearly defined 
although it is evident that the effects are not large for the 
range of Mach numbers investigated. It is expected that the 
damping would increase significantly for the attached flow 
case at 1.6 > M > 1.0. Muhlstein (ref. 18) has shown total- 
sys tern -damping measurements in attached flow at Mach 
numbers between 1. 1 and 1. 4 that decrease a full decade 
between M = 1. 1 and M = 1.4. His results at M = 1.4, 
although for a much thinner boundary layer relative to the 
panel length, are consistent with the results in Fig. 19. 

3.2.2 Displacement Response 

Illustrative measurements and computations of 
the power-spectral densities of displacement of the 
0. 001 18 -m thick panel are shown in Fig. 20 for M^ = 2.5. Measurements are shown for the probe locations indicated 
in the inset panel sketches. Corresponding computations for the same locations are shown only for the peak displace- 
ments at measured resonance frequencies. The computer program computes the complete spectrum, but this has been 
omitted to simplify the figures. The results show similar trends between attached and separated flows in the response 
of most modes as a function of frequency. The first mode predominated the displacement, as would be expected. The 
relative displacements between attached and separated flows are in the same proportion as the input spectra (Fig. 17), 
Thus the spatial correlation differences between attached and separated flows have only a subtle influence on the 
response characteristics. It was originally planned to present longitudinal and lateral acceptances as a function of 
frequency to illustrate spatial correlation effects, but these results were not available in time to be included in the paper. 

The comparison between measured and computed displacement spectral densities at resonance generally 
shows similar trends with frequency with the computed displacements from about two to ten times higher than the meas- 
ured displacements. Although the comparative results are considered good and within state of art for random excita- 
tion and response phenomena, the anxious and optimistic investigator always hopes for better agreement. The first 
suspicions are usually directed to searching for computational errors. In this case the computer program and measured 
spectra have been repeatedly checked. As an example the comparable measured spectra have been obtained by three 
different analysis systems (2 digital and 1 analog). It is believed that the major differences in computed and measured 
results are due to the assumptions in the analytical method that the flow field is uniform and homogeneous. An obvious 
major source of panel response asymmetry is the longitudinal pressure variation on the panel in separated flow 
(Fig, 4(b)). A slight lateral variation in static pressure was also found on the wind-tunnel wall in the region of the panel 
for the attached flow case. These effects are somewhat accounted for, but not entirely, by computing the response spectra 
at the measured natural frequencies. The effects of the strain-g ate s mounted on the panel are also not accounted for. 

3.2.3 Strain Response 

Typical measurements and computations of strain response power spectral densities are shown in Fig. 21. 

The most significant feature of the strains is that the first mode is not the dominant mode as shown in the case for the 
displacement response. This result is expected since the strains are proportional to the local panel curvature asso- 
ciated with the response mode. In the displacement response the differences between the attached and separated flow 
are consistent with the difference in the excitation spectra. There is also about the same order of differences between 
computations and measurements of strain. A few characteristics of the spectra indicate panel distortion and asymmetry 
of modes as discussed in 3.2.2. An example is the measurement of a significant strain response for the 1-2 mode at 
the panel center point which should be on a node line. 

3. 3 Effect of Flow Field on Displacement Response 

Comparisons of the displacement responses of two panels (d = 0 . 00 11 S m and d = 0 , 00235 m), due to the 
excitations of attached and separated and mixed flows at = 2.5 are shown in Fig. 22 , The results show that the 
response amplitudes were approximately the same for both the separated flow and mixed flow cases even though the 
excitation spectra are significantly different in the region of the shock wave (Fig. 5). The increase in mixed flow 
spectra at f < 100 Hz is consistent with the higher excitation spectra at these low frequencies. It is interesting that 
there was apparently no higher degree of coupling of the shock oscillations in mixed flow with the panel response than 
occurred in the fully separated flow case. For the mixed flow case, the shock wave was positioned near the center of 
the panels. The extent of the shock wave oscillations relative to the panel length can be seen in the longitudinal distri- 
bution of the fluctuating pressures shown in Fig. 4. 
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Figure 20. Power spectral densities of displacement response. 

The ratio between the attached flow and separated flow response spectra are approximately the same for 
both panel thicknesses, as would be expected; except that an unexplained broadband mode of vibration occurred at 
f « 100 Hz on the d = 0. 00235 m panel. Significant differences in resonant frequencies can be noted for each of the 
flow cases. The frequencies also do not scale directly with thickness or with the same scale factor for all modes for 
the two panel thicknesses. These characteristics can be due to stiffness changes resulting from un symmetrical load- 
ing or midplane stresses due to thermal effects, or due to changes in excitation; however, the distortion of the mode 
shapes by the unsymmetrical loading on the panels as discussed in 3.2.2 is considered to be the most probable cause. 


4. CONCLUDING REMARKS 

A large amount of experimental data have been studied to determine the characteristics of sui’face -pres sun 
fluctuations underlying supersonic attached and separated turbulent boundary layers (M w = 1. 6 to 3. 5). The most effec- 
tive parameters for making the characteristics dimensionless have been established. Empirical formulae have been 
derived to represent the nonsteady pressure fields in each of the flow regions. A method of analysis of the amplitude 
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Figure 21, Power spectral densities of stress response. 


and strain response of flat clamped -edge panels has also been described, and comparisons have been made between 
response measurements and computations for panels underlying both the attached and separated flows. 

The results of the pressure -fluctuation studies have shown that the turbulence generated within attached or 
separated flow is relatively homogeneous although some nonuniformities of the fluctuating pressures have been identi- 
fied. The statistical characteristics that describe the surface fluctuating pressures can be effectively scaled by the 
reduced frequency parameters f6/U or f6/U m . An attenuation- coefficient function a 6 representing the decaying 
exponential of the moduli of the cross spectra has been shown to be an appropriate function to describe the spatial 
correlation of the unsteady pressure fields. Observed irregularities in the shapes of power spectra and attenuation 
coefficients and convection velocities as a function of frequency can be attributed to changes in the predominant sources 
of the fluctuating pressures in different frequency zones. 

The results of the panel response studies indicate that the method of analysis predicts the relative levels 
of response of most modes for either attached or separated flow, but generally, over predicts the mean-square ampli- 
tudes or strains. It is believed that the differences between computations and measurements result primarily from 
the assumption in the analysis that the flow fields are uniform and homogeneous, whereas, in the experiments signifi- 
cant nonuniform steady-state loading occurred. The differences in panel response between attached and separated 
flows are consistent with the approximate three decade change in. the excitation spectra. The response of panels to the 
excitation of mixed flow with an oscillating shock wave near the center of the panels was not significantly different than 
the response to separated flow even though the excitation in the region of a shock wave has higher amplitude and differ- 
ent frequency content. Indications are that the shock oscillations did not couple to any degree with the panel motions. 
The damping of panel -airstream system was relatively invariant with Mach number at constant total pressure between 
M ffl =1,6 and = 3. 5. The damping of the first resonant mode was about 25% higher for separated flow than for 
attached flow, but otherwise the damping of all modes was about the same for either flow condition. 
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Figure 22. Comparison of panel response due to excitations of attached and 
separated boundary layer and oscillating shock waves. 
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